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We present the first 2+1 flavour spectrum measurements of glueball states using high statistics 
^-^ ' simulations with improved staggered fermions. We find a spectrum consistent with quenched mea- 

surements of scalar, pseudoscalar and tensor glueball states. The measurements were made using 
5000 configurations at a lattice spacing of 0.123 fm and pion mass of 280 MeV and 3000 configura- 
tions at 0.092 fm with a pion mass of 360 MeV. We see some evidence of coupling to 27r states. We 
compare our results with the experimental glueball candidate spectrum as well as quenched glueball 
estimates. 
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I. INTRODUCTION 

^: 

(-H ' During the last decade, a major first-principles calculation of the hadron spectrum has been carried out based on 

I— I, the improved staggered fermion discretisation formalism [l|, y|. Comparisons with experiment, have in most cases 
shown impressive agreement, typically within a few jaercent accuracy. See for example [3| and references therein and 
\^ \ a recent review article by the MILC collaboration [4^. Noticeably absent from this impressive list are some singlet 
J^ ■ quantities such as those contributing to glueball states and to the 77/77' system. Since the latter states are stable 
f — . I under the strong interaction, they may be classified as 'gold-plated' quantities Q and so deserve urgent consideration 
■n^ . as a test of this otherwise highly successful formalism. This is made all the more pressing in view of the continuing 
04 ' controversy over the validity of the so-called 'fourth-root trick' for dealing with the spurious taste multiplicity. See 
(/-v , [5| for a comprehensive review of this. A major difficulty in simulating singlet quantities is the need to make high 
f*^ ■ quality stochastic estimates of disconnected correlators which are typically very noisy observables. In dynamical 
f^ simulations the problem is particularly severe due to the high simulation costs. An analysis of these measurement 
'~H problems and initial results using MILC configurations [6| has recently been presented [3]. We have now completed 

significantly higher statistics simulations at two lattice spacings primarily in order to study the rj/rj' . (These results 
will be presented elsewhere.) As a by-product, we have also made measurements of scalar meson operators and 
glueball operators. The latter form the basis of the present paper. 

In the quenched approximation, comprehensive studies of the glueball spectrum have been available for some time [8|- 
llOJ I . In comparison, studies with dynamical quarks are still at a preliminary stage |lll4l4| . This is mainly because glue 
correlators are noisy and the masses relatively high so that signal to noise is hard to control without large statistics. 
The quenched studies typically used thousands of configurations rather than the 'few hundred' configurations usually 
available at a given lattice spacing and quark mass. We therefore took the opportunity to make use of several thousand 
configurations generated as part of the main sin glet simulation program. Preliminary results were presented in jl4l] . 
Crede and Meyer [l^ and Klempt and Zaitsev [1^1 have recently reviewed the past and future experimental searches 
for glueballs. It is particularly timely to study the glueball spectrum using unquenched QCD, because there are new 
experiments starting, or starting soon, that are looking for glueballs. The BES-III experiment IT] has begun running 
and can search for glueballs via radiative J/ip decays. It will study light mesons with ^, 2++, and 0++ quantum 
numbers, where glueball degrees of freedom could contribute. 

It will be particularly interesting to see if the existence of /j(2220) as a state is confirmed by BES-III [17]. The 
/j(2220) is a candidate 2++ state that some speculate is a glueball state because its mass is close to the mass of the 2++ 
glueball, (2390(30) (120) MeV, in the quenched estimates of Chen et al. [1^. The 120 MeV error is from different ways 
of setting the lattice spacing in quenched QCD, parametrised by uncertainty in rg. Recently the HPQCD determined 
rg to 0.8% accuracy from unquenched lattice QCD [1^. Using HPQCD's central value of l/ro = 423 MeV, changes 



the mass of the 2++ gluebah to 2470(30) MeV. This shift is within the errors of the quenched estimate [l3|, but 
ilhistrates the need for unquenched calculations. 

The PANDA experiment at FAIR will look for glueballs in the range 2.2 to 5.5 GeV [Uli^ after 2015. 

The general plan of the paper is as follows. In section|lI]we describe the simulations, the ensembles and the methods 
used for extracting spectrum information. The next section lllll contains the results from these calculations. The paper 
ends with a discussion of these results and sonic conclusions in section IIVI 



II. SIMULATION AND MEASUREMENTS 



A. Configuration ensembles 



Available resources allowed us to produce an initial two ensembles corresponding to existin g M ILC ensembles, but 
with significantly greater statistics and making use of an exact algorithm, the RHMC |2l|, l22| rather than the R 
algorithm. The simulation parameters are summarised in Table ID 



ensemble Nf 13 L^ x T 



coarse 
fine 



2+1 6.75 24-' X 64 
2+1 7.095 32^ X 64 



arrti/s 



ro/a 



a ffrnl 



0.006/0.03 3.8122(74) 0.12250(24) 
0.00775/0.031 5.059(10) 0.09230(19) 



A"cfs Ntr. 



5237 31422 
2867 17202 



TABLE I: Ensembles generated for flavour singlet studies. 



We used an improved version of the RHMC algorithm [22] which made use of such things as the 'n^'^-root trick', 
higher order integrators and different types of mass-preconditioning. Additional parameters to be tuned included 
gauge/fermion step sizes and the conjugate gradient tolerance used in different parts of the force calculation and 
acceptance tests. Details of the tuning procedures and other properties of the simulations are given in [2j]. 

The lattice spacing was determined via the static potential using standard methods [2J] . When required, the value 
of ro used to convert to physical scales was 0.467 fm. ^ 

The 'coarse' ensemble summarised in Table U represents some 8 times the statistics of the corresponding 20"^ x 64 
MILC configuration studied in [7| while the fine ensemble represents an increase of around 5 on a typical MILC 
ensemble. We have studied the autocorrelation time for simple operators. For technical reasons the first simulations 
(coarse ensemble) which were conducted on the UKQCD QCDOC machine [2^ were interrupted at various times 
resulting in non-contiguous RHMC trajectories. The ensemble was produced in two separate streams. The fine lattice 
simulations however resulted in an uninterrupted Markov chain, so allowing investigation of some autocorrelations. 
Table |ll] shows the integrated autocorrelation times of some relevant observables which can be meaningfully defined 
on a single configuration. The effective masses (defined below) are those using both local {L) and fuzzed (F) sources 





LL 


LF 


FF 


<« 


2.28(33) 


2.29(32) 


1.99(21) 




j = 


i = 1 


i = 2 


A + + 

vt' (0) 


0.65(6) 


0.89(21) 


1.12(9) 


E.n^^^i) 


1.13(15) 


1.61(22) 


2.59(26) 



TABLE II: Integrated autocorrelation times (in units of 6 trajectories) for the effective mass of the pion (averaged over euclidean 
times 8 — 10) and for the plaquette based glue operators-defined below for a fixed time plane {t — 0) and also averaged over all 
time planes. 

and sinks for the 75 x 75 pion operator. The glue operators are defined in the next section. The labels i — 0, 1,2 
refer to the number of Teper blocking levels as described below. As expected, the pion shows a longer autocorrelation 
than a single time plane operator. We have not detected any evidence of autocorrelation in the glueball correlators 
and corresponding effective masses which, as discussed below, are very noisy. We have checked, by binning in selected 
cases, that the statistical error estimates made using configurations separated by 6 trajectories are reasonable. 



A recent accurate determination by the HPQCD collaboration using MILC ensembles gave ro = 0.4661(38) [IS 



The coarse ensemble has been successfully incorporated by the MILC collaboration into an analysis of the strange 
quark content of the nucleon [2g] . 



B. Glueball measurement methods 



We used standard glue operators V^^ built from spatial plaquettes for the Oh irrep Ai~^ (coupling to 0+"*" in the 
continuum). Prior to measurement we use APE smearing [27] (twice with smearing parameter c = 2.5) and then Teper 
blocking [28[ performed n times where n = 0, 1, 2, . . . A^ — 1. Using all N blocking levels provides an A^-dimensional 
basis for the variational measurement techniques as described below. The APE smearing smooths out some of the 
ultra-violet noise while the different Teper blocking levels provide a basis of operators with a range of physical extents 
allowing different couplings to the ground state and excited states in a given channel. This is important when using 
variationally motivated techniques as outlined below. We also use these plaquette operators V^ for the E^'^ and 
T^^ irreps (7^) which couple to 2++ states in the continuum. For example, for E^^ we use both the operators in ([1]) 
(standard) and ([2|) (alternative). 



Vf^\t) = tv 



E 



(A,3;y(^iO ~ Pi,yz{^ii)) 



(1) 



and 



rf^ it) = tr ^ {Pl^y [x, t) + Ply^ - 2Pl,^ {x, t)) . 



(2) 



For pseudoscalar operators (A^ "*") we require something non-planar with some 'handedness' [23]. In practice, see 
([3|), we use left (L) and right (i?) handed versions of the loop shown in Fig[TJ 



(3) 



X 

- {Hlj.,^y{t) + H'sny^{t) + Hlj^.^^{t)) 
We refer to operators T-C^ based on these as 'hand' operators in what follows. 



"X 



FIG. 1: Schematic of the (left) hand operator "Ho which can be used to study the pseudoscalar glueball. 



We calculate correlators (vacuum-subtracted where necessary) between operators with definite momentum on time 
planes separated by euclidean time t 



C^it) = {Of{Q)Of{t)) 



(4) 



Here TZ refers to the Oh irrep under consideration and «,j refer to the operator basis being used e.g. different Teper 
blocking levels. These matrices of correlators form the basis of subsequent analyses as described below. 



1. Effective masses. These are defined generically by 



arries = In 



Cjt) 
C{t + 1) 



(5) 



where C{t) is Cjj{t) for some choice of TZ and i = j- Typically, a plateau in arrioff (i) is observed at lower values 
of t when i corresponds to higher levels of blocking. Weighted averages over the plateau region can then be 
used to estimate the ground state mass in channel TZ. When using an improved gauge action containing links 
over more than one time-plane, as here, one should be alert to possible positivity violations at small euclidean 
time J3C|. However, we have not encountered any unusual behaviour in this regard. 



2. Variational methods. In the basic variational method [31|-|33j one finds eigenvalues of the N x N matrix 

N 

Mjfit,to)^Y.iC-^(*o)]7>}C-^,{t) (6) 

fc=i 

where to is some initial euclidean time at which the correlator matrix is sufficiently well-determined to be 
invertible. It is straightforward to show (see for example the factorising fit parametrisation given below) that 
the eigenvalues X°'{to,t) of M^(io,0 ^^^ related to the transfer matrix and 

A"(to,i)=e-™"(*-*«) (a = 0,1,... TV -1). (7) 

One can then either 

(a) use directly the masses am" obtained from ([7]) or 

(b) form an effective mass, as in (O above, from the ground state {a = 0) projection of C{t) . 

For convenience, we refer to masses estimated via (a) as 'eigenvalue masses' and those estimated via (b) as 
'variational effective masses'. 

3. Factorising fits. Using the usual intermediate state arguments, the correlator matrix Q can be expressed as an 
infinite sum of exponential contributions and fitted in truncated form so as to extract mass estimates: 

Af-l 

Q,(i)=^cfc^"e-""*. (8) 

Here we have suppressed the irrep label TZ and restricted the sum to the M lowest-lying states. Note that by 
choosing M — N, the dimension of the operator basis, one can simply recover the variational formula given 
above. The factorised coefficient cf gives the overlap of state a with operator i. 

Note that in this study we have included only glue-based operator correlations in our factorising fits and correspond- 
ing variational analyses. In section iHll we will comment further on the prospects for studies of two meson operators 
and decay studies in general. 

III. RESULTS 
A. Results for the scalar glueball — coarse ensemble 

We first give some sample results of the above methods applied to the coarse lattices. Fig.[3]shows the basic effective 
masses of the momentum scalar state for each of the blocking levels to 3. Despite the significantly larger statistics, 
compared with previous dynamical simulations, it is clear that the ratio of signal to noise is still a significant problem. 
Weighted averages of the 'plateau' values are shown for blocking levels 2 and 3. Here the weighting is inversely 
proportional to the statistical error. As with all other quantities in this study, overall statistical errors were estimated 
via bootstrap. 

Fig. [3] shows the variational effective masses, as defined in section Hi Bl (method 2b) deduced from a 3 x 3 correlator 
matrix (blocking levels to 2) 

Similar effective mass estimates were made using momentum 1 operators. The latter require no vacuum subtraction. 

In Fig. |4] we show the lowest 3 masses (where obtainable) extracted from the variational matrix (method 2a). 
The momentum 1 masses have been obtained from energies using the nai've lattice dispersion relation. There is a 
reasonable degree of consistency with respect to the choice of t/to except at large t where there is a suggestion that 
the lowest state projected out corresponds to a tttt system. Given that for these light quark masses, the two n decay 
threshold is open, this is not unexpected. 
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FIG. 2: The efTective masses for the diagonal entries (0, 0), (1, 1) etc. of a 4 x 4 matrix of correlators formed using the |p| — 
standard scalar operators with blocking levels to 3 measured on the coarse ensemble. Weighted averages for blocking levels 
2 and 3 are shown, computed from f = 2 — 5. 



Similar studies were made using the alternative hand operators defined in section |TTbJ The results were statistically 
consistent with those obtained from the standard operators. Fig. [5] shows a comparison of the lowest three states 
using standard and alternative operators. 

Finally in this subsection, we present results of the third method: multichannel factorising fits. We have carried 
out 3x3 and 4x4 fits with 2 and 3 exponentials. We have studied both correlated and uncorrelated fits. We also 
investigated stability with respect to the fitted t range tmin l£t ^ ^max. For example. Fig. |B] shows the variation with 
^min of the lowest two states resulting from two and three exponential fits to a 4 x 4 matrix of correlators (blocking 
levels — 3). These fits are fully correlated and the corresponding fit details are listed in Table Hill 

The extent of consistency between the various methods can be judged from Fig.[71 The solid line, with dashed errors, 
represents a global average of these determinations. This average was obtained using the methods adopted by the 
Particle Data Group [3J| when combining results from independent determinations. The errors used in the weighting 
are combined systematic and statistical errors. Of course, the statistical errors in our case are not independent. 
For each type of determination contributing to Fig. [71 represented by a central value and error bar, we have used 
the same weighted averaging procedure. The systematic errors took into account the variation due to fit range. In 
fact, the determinations contributing to the global average are a selected sub-group of all the various estimates and 
fits described above [2J|. We have excluded fits which had unacceptably high (or low) x^/d.o.f. and/or which had 
very large errors, and also those which were found to have high sensitivity to fit ranges and/or parameter choices. 
Similar considerations were applied to the determination of the first excited state in the scalar sector, also shown in 
Fig. [71 We have only quoted fitted mass values for states where at least one higher state has been included in the 
determination. The global averages for the lowest two states in the scalar channel are shown in Table IIVI along with 
the corresponding values for the fine lattices (described in the next subsection). We will return later to the physical 
interpretation of these results. 
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FIG. 3: Effective masses computed using tiie variational correlator for different choices of t/to- 
for the 2/1 and 3/2 projections (coarse ensemble, standard scalar glue operators). 



Weighted averages are shown 



iVexp 


\P\ 


(^niin 


(^max 


arno 


ami 


am2 


xVd.o.f. 


2 





1 

2 
3 


6 
6 
6 


1.0336(29) 
0.887(14) 
0.62(13) 


1.84(11) 
1.332(47) 
1.22(10) 


— 


3.176 
0.635 
0.349 


2 


1 


1 
2 

3 


6 
6 
6 


1.0770(60) 
1.034(36) 
1.19(15) 


1.70(14) 
1.42(22) 
3.30[-] 


— 


0.832 
0.301 
0.189 


3 





1 
2 

3 


6 
6 
6 


0.9409(43) 
0.32(15) 
0.306(92) 


1.626(37) 
0.9495(57) 
1.016(17) 


18.04[-] 
1.77hl 
7.79[-l 


1.309 
0.607 
0.327 


3 


1 


1 
2 


6 
6 


1.0682(60) 
1.046(31) 


1.623(75) 
1.41(20) 


21.26H 
3.53[-] 


0.329 
0.173 



TABLE III: Fitted mass parameters for two and three-exponential factorising fits to a 4 x 4 matrix of correlators using blocking 
levels — 3 (coarse ensemble). Where errors are quoted as [— ] this indicates that the gradient in that direction of parameter 
space was undetermined. 



ensemble 


a [fm] 


am 


am. 


coarse 
fine 


0.12250(24) 
0.09230(19) 


1.0468(75) 
0.8332(59) 


1.875(87) 
1.368(17) 



TABLE IV: Global average values for scalar masses on coarse and fine ensembles 



B. Results for the scalar glueball — fine ensemble 



Fig. [8] (corresponding to Fig. |3] for the coarse ensemble) shows variational effective scalar masses and weighted 
averages on the fine ensemble. In Fig. ^ we show the lowest 3 masses extracted from the variational matrix as 
described in method 2a. This is to be compared with Fig. 2] for the coarse ensemble, but in this case we take the 




FIG. 4: The masses extracted from the variational eigenvalues for different f/fo using blocking levels 0, 2, 3 for momentum 
and 1,2,3 for momentum 1. The energies corresponding to a mr state with momentum and 1 are drawn for comparison 
(coarse ensemble, standard scalar glue operators). 



opportunity to display results for both standard and alternative (hand) operators. Both operator bases gave consistent 
results. 

As examples of factorising fits on the fine ensemble, we present the imax dependence at fixed imin hi Fig. 1101 These 
are 4x4 fits using two and three exponentials. We see excellent consistency for the groundstates between the M = 2 
(both \p\ — and 1) and the M — 3, \p\ = 1 fits. The corresponding fit quality can be judged from Tables IVl and I VII 

As with the coarse ensemble, we have selected the most reliable determinations for each of the described methods 
and assigned systematic errors to cover sensitivity to parameter and fit range. As an example, Fig. 1111 shows sample 
variational results (efTective mass and eigenvalue methods) for momentum zero and one. The relevant time slices t/to 
are shown on the horizontal axis along with the t range used for averaging the effective mass. 

The final values selected for all methods are shown in Fig. [T^] along with the resulting global average using the same 
weighting procedures described earlier in section IIII Al 



C. Scalar decay and mixing 



Having observed signals in our fits which appear close to the 2tt threshold, we have attempted to check the strength 
of mixing between the glueball and tttt operators. Early attempts to study the mixing directly have been made 35 1 
but were performed at pion masses such that 2mT^ was just below ruQ. One might naively expect that glueball decay 
would be flavour blind, decaying into tttt, KK and rjrj final states with equal rates. Of course, since we are now 
observing the state (or states) to which the gluonic operators couple most strongly, then mixing effects are likely to 
violate flavour blindness. If the 'glueball' mixes strongly with a mainly uu + dd state then the tttt decay channel is 
preferred over the OZI suppressed KK channel, and vice versa if the coupling is strongest to an ss state. There are 
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FIG. 5: The masses extracted from variational eigenvalues for different t/to (coarse ensemble). These were made using blocking 
levels 0, 1,2 for standard plaquette (filled circles) and alternative hand (empty circles) operators. 



chiral suppression arguments [36| that suggest that KKbar could dominate the two meson decay channel in a glueball 
state in this mass range. 

One can in principle study decay matrix elements directly on the lattice following the procedure of Lellouch and 
Liischer [37| . However this requires very accurate determinations on different lattice sizes and is a difficult technique 
to apply even for states which are much less subject to noise than those in the glueball sector. Initial attempts to 
study resonant states with open decay channels have been made - for example the PACS-CS collaboration [38] has 
studied the isovector P-wave phase shift in tttt scattering in order to estimate the width of the p meson. Additional 
complications and extra diagrams introduced in the case of staggered fermions have been detailed for example in [3^ 
where 1 = 2 phase shifts were studied. In the isoscalar channel there are of course the further complications arising 
from disconnected contributions. The possible types of correlators required for a mixing study of glueball, /o and tttt 
are indicated in Fig. [131 Unquenched studies of mixing in the isoscalar 0++ system were carried out by Hart et al. [ij] 
using two flavours of Wilson fermions. The ratio of diagram (i) to to diagram (h) (see Fig. [T5)) was estimated using all 
to all techniques with 100 configurations. Unfortunatly the errors were too large to provide useful information within 
the statistics of that study. 

It has become clear that techniques for studying lattice decays will best be developed first in simpler and less noisy 
systems than isoscalar channels such as the present one. 

So, for a preliminary study with improved staggered fermions, we choose instead to follow a simplified version of 
the method of [4Q| [ and form the mixing matrix 
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{vt\mf\t)) {CAQ)rf\t)) 



(9) 



where P^ ^ (0) and P^ ^ (t) are the \p\ = standard scalar glueball operators at timeslices and t respectively, 
and C7r(0) and CTr(i) are the single pion correlators on timeslice and t respectively. The C7r(0) correlator allows 
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FIG. 6: Results from factorising fits applied to the 4x4 basis of standard scalar glueball operators with blocking levels to 3. 
Both two and three-exponential fits are presented for momentum and 1 correlators, with fmax ~ 6 and fmin allowed to vary 
(coarse ensemble). 



us to study a tttt state localised in time and, because of the way we have computed the connected correlator, we are 
restricted to using CtAO) which reduces our statistics for the 27r operators by a factor of n* = 64. 
We form the ratio 1401 



XGirit) 



Ma^it) 



^ MGG{t)M^^{t) 



(10) 



where the Mab{t) are the elements of ^ on timeslice t. This gives a measure of the off-diagonal mixing matrix 
elements, normalised by the diagonal entries. We remind the reader that this captures only some of the possible 
contractions contributing in a full treatment of the mixing as indicated in Fig. 1131 In this study, the factorising fits 
involved diagrams of type (a) only and the estimated strength of mixing to tttt involved an approximation to (c) . 
We have computed these ratios for the coarse and fine lattices, using the glueball operator with three levels of Teper 

blocking (i.e. V^ ^ ) which we found shows similar behaviour to the lower blocking levels but with less noise. We 
used the pion correlators with both local source and local sink, and fuzzed source and fuzzed sink. Our results are 
presented in Fig. [141 

For the coarse ratios we see that at small t they appear consistent with zero, turning negative for i ^ 5, although 
with large error. We note that the fine mixing ratio shows a similar downturn for i '^ 3 — 4. Whilst these ratios 
give only a guide to the mixing of the glueball operators and a tttt state, the small size of xgtt at low t indicates that, 
provided we choose t^-m small enough, we should obtain a good overlap with the 'stable' glueball and conversely, 
that by choosing t^in large, we may obtain a significant overlap with the tttt state. This is consistent with the earlier 
observations of the effective masses and factorising fits at small and large t. 
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FIG. 7: The overall average for the scalar glueball groundstate and first excited state masses on the coarse ensemble, computed 
as described in the text. 



D. Comparison with other scalar determinations 



This study of scalar states in the glueball sector is characterised by relatively high statistics in comparison with 
previous studies using dynamical fermions. It also uses a vacuum with 2 + 1 flavours. However, it is limited in that 
it makes use of only two lattice spacings (roughly 0.12 and 0.09 fm) and a single light quark mass for each ensemble. 
The pion mass in lattice units is 0.1740(6) (coarse ensemble) and 0.1672(14) (fine). These correspond to around 280 
MeV and 360 MeV respectively so the quark masses are not particularly light. The strange quark mass is close to the 
physical one. Any attempt to extract a continuum limit (without a chiral limit of course) needs to be treated with 
due caution. In Table IVTIl we summarise the scalar sector results in lattice units and in dimensionless form using the 
Sommer rg parameter from Table HI 

In Fig. [15] we present our results along with those from a recent UKQCD study of scalar glueballs using staggered 
fermions (MILC configurations) 41], a UKQCD study using 0{a) non-perturbatively improved Wilson fermions 
{Nf = 2) [111] and the continuum limit result from a quenched anisotropic study of the glueball spectrum ^9]]. We 
perform a simple continuum extrapolation using the form 



rom(a) ^ romc + b{rQ/aY 



(11) 



where r^mQ and b are the parameters to be determined. In [9], this was found to work well for all glueball states 
except for the A^^ where this was thought to be due to strong lattice spacing dependence caused by proximity to 
the non-physical phase transition in the Wilson fundamental-adjoint plane for the anisotropic lattice action used. 
However, the alternative form employed in [9| has four free parameters which we are unable to use with just two data 
points. As noted above, it is ambitious to attempt any kind of continuum extrapolation using two data points only and 
we emphasise that the fit with pTjl has been performed as a guide only. Nevertheless, we that note good consistency 
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0.8446(13) 


18.93[-] 
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0.8480(11) 


1.54H 
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1.933 



TABLE V: Fitted mass parameters for two-exponential fully correlated factorising fits to a 4 x 4 matrix of correlators using a 
blocking levels 0, 1,2,3 with \p\ = and 1 for varying tmax with imin fixed at t = 1 (fine ensemble). Where errors are quoted 
as [— ] this indicates that the gradient in that direction of parameter space could not be determined. 



between our continuum value for romc so obtained, 4.32(6), and the quenched value romc = 4.21(11)(4) 9]. Our 
extrapolated value corresponds to a physical value of 1.83 GeV approximately . 

Our masses seem to show rather weak dependence on a — certainly weaker than that observed in |ll| and of a 
similar strength to that observed in [9\ . 

In Fig. [in] we present our results plotted against the pion mass, shown with the same comparisons from the literature 
as in Fig. [15] One might tentatively claim that the glueball mass increases as one decreases the pion mass, hinting 
at some underlying mixing dynamics. However the UKQCD measurements on the coarse MILC Asqtad ensembles 
are rather spread out and if^as is suspected, the 0{a) improved Wilson measurements are suppressed by the phase 
structure of the action used [ll| then they should probably be discounted and there remains very little trend to study. 



12 




FIG. 9: Scalar glueball masses extracted, for the fine ensemble, from the variational eigenvalues using blocking levels 0, 1, 2 for 
different f/to. Results based on the standard and alternative (hand) operators are compared. The energies corresponding to a 
TTTT state with relative momenta \p\ — and |p| = 1 are drawn for comparison (fine ensemble). 



E. Tensor and pseudoscalar glueball results 



Similar procedures were used to obtain estimates of the ground state and first excited state in the pseudoscalar 
and tensor glueball channels. For the former we used blocked operators (jH]) based on the H operators while for the 
latter we used those based on P - see ([T]) and ^. Sample eigenvalue masses are shown in Fig. [TTland Fig. [TS] for the 
pseudoscalar and tensor states respectively. 

The global averages are given in Table IVIIII For technical reasons, we were unable to complete measurements of 
the tensor state on the fine ensemble. For comparison, the quenched continuum limit estimates of rom from [9| are 
6.33(13) and 5.85(8) for pseudoscalar and tensor glueballs respectively. So, just as for the scalar channel, there is 
little evidence of strong unquenching effects. 

We have conducted a separate study of the rj/i]' system using connected and disconnected meson operators which 
will be presented elsewhere. Just for orientation, we note that the continuum (experimental) values of rom would be 
1.32 (2.27) for the rj [r]') respectively. It is clear that there should be little mixing between the pseudoscalar glueball 
and the rj' . However, further excited states in this singlet channel may mix. 

There have been claims that the pseudoscalar ?7(1405) meson [4^,|43l contains significant mixtures of the pseudoscalar 
glueball. This requires that there is a significant shift of the mass of the pseudoscalar glueball between quenched 
and unquenched QCD, which we don't find in this calculation. However, note the qualifications accompanying our 
conclusions below. 
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and tmax allowed to vary. 



IV. DISCUSSION AND CONCLUSIONS 



We have presented the first high statistics study of the low-lying glueball sector in 2 + 1 flavour QCD. We analysed 
some 5000 configurations of a 24^ x 64 lattice with spacing 0.12 fm and 3000 configurations of a 32^ x 64 lattice with 
spacing 0.09 fm using improved staggered (Asqtad) dynamical fermions. In contrast with earlier glueball analyses 
using improved Wilson fermions, we find no evidence of strong unquenching effects. For the 0^"*", ^ and 2"*"+ states 
we obtain mass estimates quite close to the continuum limit masses obtained in the quenched approximation [9|. In 
Fig. [m we show a compilation of 0++, ^ and 2++ glueball candidate states [3J) together with quenched continuum 
limit predictions [9| and our unquenched results at fixed lattice spacing and fixed quark mass. 

Despite the larger statistics available to us, we still had difficulty obtaining unambiguous estimates of even the 
low-lying states. This was due to at least two main factors - the inherently noisy correlators and the presence of open 
decay channels, particularly tttt. In order to reduce the noise, we used a variety of smearing and blocking techniques 
giving access to a range of basis states. We also used a number of different methods based on variational techniques: 
effective masses, transfer matrix eigenvalues and multi-channel factorising fits. The effects of the open decay channels 
were exposed at larger euclidean time via the choice of fit ranges and also via direct evaluation of a subset of mixing 
matrix elements. 

Our overall conclusions are: 

• there is little evidence of large unquenching effects on the predicted low-lying glueball spectrum; 

• accurate determination of masses requires even larger numbers of configurations in comparison with quenched 
glueball analyses; 

• a fuller account must be taken of two meson state contributions at large euclidean time; 
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TABLE VI: Fitted mass parameters for three-exponential fully correlated factorising fits to a 4 x 4 matrix of correlators using 
a blocking levels 0, 1, 2, 3 with |p| = and 1 for varying tmax with tmin fixed at t = 1 (fine ensemble). 



Result 


am(yl++) 


rom(A++) 


Coarse - Ground 
Fine - Ground 

Coarse - Excited 
Fine - Excited 


1.0468(75) 
0.8332(59) 
1.875(87) 
1.368(17) 


3.991(36) 
4.215(38) 
7.15(35) 
6.92(10) 



TABLE VII: Scalar glueball masses (ground and first-excited states) from the coarse and fine lattices converted into units of 
the Sommer parameter tq. 

• we encountered no problems that could be identified as resulting from the use of improved staggered fermions 
- in particular, the lattice spacing dependence was weak. 
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TABLE VIII: Pseudoscalar and tensor glueball masses (ground and first-excited states) from the coarse and fine ensembles 
converted into units of the Sommer parameter ro. Pseudoscalar measurements were made on 3506/1998 configurations of the 
coarse/fine ensembles respectively. Tensor measurements were made on 2627 configurations of the coarse ensemble only. 
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FIG. 12: The overall average for the scalar glueball groundstate mass on the fine ensemble, evaluated as described in the text. 
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FIG. 13: Schematic mixing diagrams for the full glue, meson, and two meson mixing problem in the / = scalar channel. 



The dashed rectangles represent the glue operators V 
horizontally. 



and the solid lines represent quark propagators. Euclidean time runs 
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FIG. 14: The measure of giueball-TTTr mixing defined in (|10p for, clockwise from top-left: the local-local (LL) pion correlators 
with the three-times Teper blocked glueball operators on the coarse ensemble; the fuzzed-fuzzed (FF) pion correlators with the 
three-times blocked glueball operators on the coarse ensemble; the fuzzed-fuzzed (FF) pion correlators with the three-times 
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FIG. 15: Our measurements of the scalar glueball mass shown with previous quenched (Morningstar and Peardon [y]) and dy- 
namical (UKQCD on MILC Asqtad [4l|] and UKQCD Clover [ll|]) determinations with the continuum extrapolation performed 
as in [^. 
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FIG. 16: Our measurements of the scalar glueball mass shown with previous dynamical (UKQCD on MILC Asqtad [4l|] and 
UKQCD Clover [TlJ) determinations. The tttt threshold is shown (dash-dotted line) for convenience. 
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FIG. 17: The masses extracted from the groundstate variational eigenvalues using different t/to projections performed on a 

A — '' A — ^ A ^ 

3x3 matrix of correlators formed from the basis of momentum zero pseudoscalar glueball operators {"Hg ^ , Tii ^ , ^-2 ^ } 
on the fine ensemble. 
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FIG. 18: The masses extracted from the groundstate variational eigenvalues for different t/to projections performed on 4 x 4 



matrices of correlators formed using bases of momentum zero tensor glueball operators P^i 
blocking levels 0, 1, 2 and 3 in each case (coarse ensemble). 
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FIG. 19: Comparison of our unquenched glue-sector predictions with experimental states listed in the Particle Data Group 
listings [3J and with the quenched continuum limit predictions of Morningstar and Peardon [y]. The mass range expected to 
be investigated in a forthcoming experiment 19i] is also indicated 



